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ABSTRACT 

The recently detected linear polarization in the optical lightcurve of GRB 990510 re- 
newed the interest on how polarization can be produced in gamma-ray burst fireballs. 
Here we present a model based on the assumption that we are seeing a coUimated fire- 
ball, observed slightly off-axis. This introduces some degree of anisotropy, and makes 
it possible to observe a linearly polarized flux even if the magnetic field is completely 
tangled in the plane orthogonal to the line of sight. We construct the lightcurve of the 
polarization flux, showing that it is always characterized by two maxima, with the po- 
larization position angle changing by 90° between the first and the second maximum. 
The very same geometry here assumed implies that the total flux initially decays in 
time as a power law, but gradually steepens as the bulk Lorcntz factor of the fireball 
decreases. 
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1 INTRODUCTION 

It is now widely believed that the afterglow emission of 
,: gamma-ray bursts is due to the deceleration of the rela- 
tivistic fireball in the circum-burst matter (for reviews see 
Piran 1999; Meszaros 1999). This produces a shock that ac- 
. celerates electrons to random relativistic energies and prob- 
' ably enhances the magnetic field, leading to the production 
of synchrotron emission. If the magnetic field is completely 
tangled over the entire emission seen by the observer, the 
resulting synchrotron emission is unpolarized. On the other 
hand a very high degree of linear polarization can be ex- 
pected if a fraction of the magnetic field is well ordered, 
reaching 60-70% in the case of a completely ordered field. 
Polarization values in the optical band in the range 3-30% 
have been indeed observed in cosmic sources, like BL Lac 
objects and High Polarization Quasars (see e.g. Angel & 
Stockman 1980; Impey & Tapia 1990), whose radiation is 
believed to be produced by the synchrotron process. One 
therefore expects that also in gamma-ray burst afterglows 
the emission is polarized, and attempts were made to mea- 
sure it. After the upper limit (2.3%) found for GRB 990123 
(Hjorth et al. 1999), Covino et al. (1999) detected linear 
polarization in the afterglow of GRB 990510, at the small 
but significant level of 1.7±0.2%. This detection was then 
confirmed by Wijers et al. (1999) who detected similar po- 
larization values two hours and one day later. 

On the theoretical side, Gruzinov & Waxman (1999, 



hereafter GW99) and Gruzinov (1999) predict values around 
10%, significantly larger than observed. This estimate is 
based on the assumption that the overall emission reach- 
ing the observer is produced in a finite number A'^ ~ 50 of 
regions causally disconnected, each of which is embedded 
in a completely ordered magnetic field. The predicted total 
polarization level is 60%/viV, equal to ~10% for A'' ~ 50. 
GW99 discuss how the coherence length of the magnetic 
field generated at the external shock front of a GRB fireball 
grows with time. If, however, the magnetic field is generated 
at the coUisionless shock front, which is extremely local, it 
is not clear why the magnetic field embedded in the newly 
swept matter should be linked to the field in the regions 
behind the shock. 

An alternative magnetic field generation process (and 
hence geometry) has been discussed by Medvedev & Loeb 
(1999, ML99 hereafter), who consider a magnetic field com- 
pletely tangled in the plane of the shock front, but with a 
high degree of coherence in the orthogonal direction. In the 
case of a spherical fireball this geometry produces no po- 
larization unless a part of the fireball emission is amplified 
and part is obscured, as is the case of interstellar scintilla- 
tion. In this case, however, the resulting polarization can be 
much better observed at radio wavelengths and should show 
a rapid and erratic change of the position angle. 

We here propose an alternative model, in which the 
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magnetic field geometry is analogous to that of ML99r]but 
in a fireball that is collimated in a cone and observed slightly 
off-axis. In this case the circular symmetry is broken and 
net polarization can be observed (see e.g. Hjorth et al. 1999, 
Covino et al. 1999, Wijers et al. 1999). Evidences for beam- 
ing of the fireball of GRB 990510 from the anomalous decay 
of the optical lightcurve has been discussed in many recent 
papers (Harrison et al. 1999, Israel et al. 1999, Stanek et al. 
1999b). 

The key assumption of our model is that the fireball 
is collimated in a cone, observed slightly off-axis. The key 
result we obtain is the polarization lightcurve, its connection 
with the flux behavior and a characteristic change of 90° in 
the polarization angle, making the model very testable. 



2 POLARIZATION LIGHTCURVE 

2.1 Magnetic field configuration 

Assume a slab of magnetized plasma, in which the config- 
uration of the magnetic field is completely tangled if the 
slab is observed face on, while it has some some degree of 
alignment if the slab is observed edge on. Such a field can 
be produced by compression in one direction of a volume 
of 3D tangled magnetic field (Laing 1980, hereafter L80) or 
by Weibel instability (ML99) . If the slab is observed edge- 
on, the radiation is therefore polarized at a level, Po, which 
depends on the degree of order of the field in the plane. 
At the angle 6 from the normal of the slab, the degree of 
polarization can be expressed by, following L80: 



Relativistic Fireball 



p{e) 



\e) 



l-Hcos2(6i) 



(1) 



If the emitting slab moves in the direction normal to its plane 
with a bulk Lorentz factor F, we have to take into account 
the relativistic aberration of photons. This effect causes pho- 
tons emitted at 6' = 7r/2 in the (primed) comoving frame 
K' to be observed at 61 ~ l/F (see also ML99). 



2.2 Polarization of beamed fireballs 

We assume that in gamma-ray burst fireballs the emitting 
region is a slab expanding radially and relativistically, com- 
pressed along the direction of motion. We assume also that 
the fireball is collimated into a cone of semi-aperture angle 
6c, and that the line of sight makes an angle 9o with the jet 
axis (upper panel of Fig. m. As long as F > l/(Sc — do), the 
observer receives photons from a circle of semi-aperture an- 
gle 1/F around 6o (i.e. within the grey shaded area of Fig. 0). 
Consider the edge of this circle: radiation coming from each 
sector is highly polarized, with the electric field oscillating in 
radial direction (see also ML99). As long as we observe the 
entire circle, the configuration is symmetrical, making the 
total polarization to vanish. However, if the observer does 
not see part of the circle, some net polarization survives in 
the observed radiation. This happens if a beamed fireball is 
observed off-axis when l/{9c + 6o) <T < l/(6'c - 6^0). 



* Note, however, that the ML99 instability is not the only process 
that can be responsible of such a geometry, see e.g. Laing 1980. 





Figure 1. Geometry of the beamed fireball (upper panel). 9o is 
the observer angle, 9^ is the cone aperture angle. The lower panel 
shows a zoom of the region around the line of sight. Note that 
the photons emitted in the comoving frame at an angle 7r/2 from 
the velocity vector are those making an angle 6 ~ l/F with the 
line of sight in the observer frame. 



The probability to observe a cone along its axis is van- 
ishingly small, since it corresponds to a small solid angle; it 
is therefore more likely to observe the collimated fireball off- 
axis. If the cone angle 9c is small, the probability p{9o/0c) 
is approximately distributed as: 

where {9o) is the average off-axis angle. 

Assume therefore Oo/Oc > (Fig. El). At the beginning 
of the afterglow, when F is large, the observer sees only a 
small fraction of the fireball (grey shaded region in Fig. El) 
and no polarization is observed. At later times, when F be- 
comes smaller than l/{dc — So), the observer will see only 
part of the circle centered in Oo'. there is then an asymmetry, 
and a corresponding net polarization flux (horizontally line 
shaded region of Fig. 0). To understand why the polariza- 
tion angle in this configuration is horizontal, consider that 
the part of the circle which is not observed would have con- 
tributed to the polarization in the vertical direction. This 
missing fraction of vertical polarization does not cancel out 
the corresponding horizontal one, which therefore survives. 

At later times, as the fireball slows down even more, a 
larger area becomes visible. When F ~ l/{6c + 0o), the dom- 
inant contribution to the flux comes from the upper regions 
of the fireball (see Fig. H), which are vertically polarized. The 
change of the position angle happens when the contributions 
from horizontal and vertical polarization are equal, result- 
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Figure 2. Front view of the beamed fireball, as in Fig. 1. C.A. 
stands for Cone Axis, while L.o.S. stands for Line of Sight. The 
grey shaded area produces unpolarized radiation, due to the com- 
plete symmetry around the line of sight. The horizontal line 
shaded region produces a horizontal component of polarization 
while the upper region (vertically line shaded) produces vertical 
polarization. 



ing in a vanishing net polarization. At still later times, when 
r ^ 1, hght aberration vanishes, the observed magnetic field 
is completely tangled and the polarization disappears. 

We therefore expect two meixima in the polarization 
lightcurve, the first for the horizontal component and the 
second for the vertical one. 

The quantitative calculation of the lightcurve of the po- 
larized fraction has been worked out as follows. We assume 
that the emission of each small volume of the fireball is 
isotropic in the comoving frame K' . We also assume that in 
this frame each small element emits the same instantaneous 
intensity. For each ring of radius p and width dp of fireball 
material around the line of sight (see Fig. H) the relativisti- 
cally enhanced monochromatic intensity has been computed 
as I{u) = 5^r{v'), where 5 = [r(l - /3cos6')]-^ is the rela- 
tivistic Doppler factor, and v = 5p' . For I' iy') oc {i'')~°' we 
then have I{v) = 5'^^°'I'{v). Fig. 3 shows the geometrical 
set up of the system. In our calculations, the distances R, 
Ri and p are considered as angles, since all of them scale 
with the distance d of the fireball from the center of explo- 
sion (see Fig. 1). In this case, p is equivalent to the angle 
d between the velocity vector of each element and the line 
of sight. Therefore the intensity observed from each element 
is a function of F and p. The total intensity is obtained by 
integrating over the entire surface of the fireball, taking into 
account that each ring is observed at a different angle, and 
then characterized by a different 5. We assume a constant 
spectral index a = 0.6 throughout our calculations. The to- 
tal intensity from the entire fireball at a given time (and 
hence at a given F) is: 




Figure 3. Sketch of the geometrical set up used to compute the 
total and the polarized flux. 



7(r,i.) = 27r I{r,u,p)dp+ 



I{r,,y,p)[2TV- 
Ri 



2Mp)]Ms) 



The first integral corresponds to the grey shaded area of 
Fig. 2, and to p < _Ri (see Fig. 3). The second integral 
corresponds to p > i?i. In this case we receive radiation 
only from a sector of the ring, between the angles tpi and 
27r — -i/ji, with il)i given by: 



tpi 



IT 

2 



2RiR -Rj- p^ 
2p{R^Ri) 



(4) 



We now compute the polarized intensity. To this aim, 
we consider again each element of the ring of radius p. After 
calculating the corresponding viewing angle in the comoving 
frame, we apply Eq. 1 to derive the intensity of the linearly 
polarized light. The position angle of the polarization pro- 
duced by each element is in the radial direction. We then 
calculate the polarization of each ring by summing the po- 
larization vectors. Finally, we integrate over p. 

It is convenient to write the polarization vector as a 
complex numbeir] PiV, v, p) = P(r, v, p)e^**'' and integrate 
it between ipi and 27r— t/ji (see Fig. 0). Here 6p is the position 
angle of the linear polarization in the observer frame. The 
polarization of a generic ring is: 

i'2-n — ^l 

P{T,v,p) = P{T,v,p) / e'^*dV = nr,^,p) sin (2,^1) (5) 



Since the result is a real number, the polarization direction 
can lie either in the plane that contains both the line of sight 
and the cone axis (we call it "vertical" polarization angle) 
or in the orthogonal plane (i.e. "horizontal" polarization an- 
gle). No intermediate values are possible. 

T This is fully equivalent to the Stokes notation of the Q and U 
parameters for linear polarization. 
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Integration over p then yields: 



P{V,u) 



n^^^)jR, 



2R~Ri 



I{V,v,p)P{V,v,p)dp I e^'^d^ 



27r — i/>i 



ipl 



2R-Ri 



. I{V,v,p)P(V,v,p)sm[2xP^[p)] dp (6) 

'^y^ ^'^1 J R^ 

For simplicity, we neglected the light travel time effects 
introduced by the overall curvature of the emitting regions 
(see Fig. hi), approximating the emitting volume with a slab. 
Since the degree of polarization of each element is divided 
by the total intensity of the same element, we expect these 
effects to be small in our case. Note that to include this 
effect requires to assign a specific relation between F and the 
observed time t, which can be different according to different 
models (due to e.g. adiabatic vs radiative evolution and/or 
gradients in the interstellar density). 

The result of the numerical integration of Eq. ra is shown 
in the lower panel of Fig. M for four different values of the 
off-axis ratio 9o/0c- For the specific cases shown in the fig- 
ure we have assumed 6c = h° , but the general properties of 
the polarization lightcurve are unaffected by the particular 
choice of 9c- All the lightcurves (except the one with the 
lower off-axis ratio, which shows almost zero polarization), 
are characterized by two maxima in the polarized fraction. 
As discussed above, at the beginning the polarization is hor- 
izontal, it reaches a first maximum and then, when the po- 
larized fraction vanishes, the angle abruptly changes by 7r/2. 
The polarized fraction raises again and then finally declines 
to zero at late times. The second peak has a value always 
larger than the first one. 

By interpolating the results of the numerical integra- 
tions (assuming a — 0.6), we can express the maximum of 
the polarization lightcurve as a function of the off-axis an- 
gle: 



p 

-t max 



0.f9P( 



■« m 



(7) 



^<^<i; l°<e.<l5° 



which is accurate within 5% in the specified 8o and 9c ranges. 
This maximum always corresponds to the second peak. 



3 LINK WITH THE TOTAL FLUX 
LIGHTCURVE 

The scenario proposed above for the polarization behav- 
ior has a strict and direct link with the behavior of the 
lightcurve of the total flux. The exact connection will de- 
pend on the specific model assumed for the deceleration of 
the fireball (i.e. homogeneous or radially distributed density 
of the interstellar medium, adiabatic or radiative regime, 
and so on), but there are general properties that can be 
incorporated in any model. For illustration, assume the sim- 
plest model of a spherical fireball expanding adiabatically 
in a homogeneous medium, predicting that the bulk Lorentz 
factor F oc i~^' *, and predicting a decay law for the flux den- 
sity J'„(t) ex f-3"/2 (Meszaros & Rees 1997), where a is the 
spectral index of the radiation spectrum, i.e. F^ oc 1/'°". This 
law assumes that the flux is proportional to the source solid 
angle: if the fireball has reached a region of size d, then the 
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Figure 4. Lightcurves of the total flux (upper panel) and of the 
polarized fraction (bottom panel) for four different choices of the 
ratio 6o/dc- The cone aperture angle 8c = 5°. Oo is the viewing 
angle, as defined in Fig. hi. The higher the ratio do/dc, the higher 
the polarized fraction due to the increase of the asymmetry of 
the geometrical setup. The actual value of the observed polar- 
ization depends linearly upon Pq (see text). For this figure we 
assumed Pq = 60%. The lightcurve of the total flux assumes a 
constant spectral index a = 0.6 for the emitted radiation. Note 
that the the highest polarization values are associated with total 
flux lightcurves steepening more gently. To calculate the values 
of the upper x-axis (t*), we assumed (t*/to) = (F/Fo)"*''' with 
to = 50 s and Fq = 100. 



accessible solid angle Q. oc (d/F) oc (iF) oc t^' . In the case 
of a fireball coUimated in a cone of constant opening angle 
(e.g. neglecting for simplicity the possible lateral spreading 
of the firebaU: Rhoads 1997, 1999; see see also Moderski, 
Sikora & Bulik 1999), we will have the above relation as 
long as 1/F < 9c — 9o- On the other hand, when 1/F becomes 
greater than 9c + 9o, all the cone front becomes visible, and 
therefore Q oc {9cdf oc {9ctr'f oc t^/^. This produces a 
steepening of the power law lightcurve decay, which now be- 
comes F„{t) oc i-3°=/2-3/4 ^ggg g^jgg Meszaros & Rees, 1999). 
At intermediate times, for which 9c — 9o < 1/F < 9c + 9o, 
the accessible solid angle increases with a law intermediate 
between i^'* and i^", producing a gradual steepening in 
the lightcurve. 

Assuming again F oc t~^'^, and calling Fi = l/{9c — 9o), 
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r2 — l/(Sc + do) and ii, t2 the corresponding times, we 
simply have 



£1 
Fa 



tl. 



(i2Al)^/« 
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(t2/tl)3/« + 1 



(8) 



In the case of the GRB 990510 afterglow, the decay laws 
have been i~°'^ at early times (Galama et al. 1999) up to 
~ half a day after the burst event, and t~^'^ after ~5 days 
(Israel et al. 1999). Therefore the ratio ta/ii is of the order 
of 10-15, implying 9o/6c between 0.4 and 0.45. However, 
the observed steepening was larger than the value (3/4) de- 
rived above. Another cause of steepening can be the lateral 
spreading of the fireball, as suggested by Rhoads (1997), 
when r ~ l/6c (but see Panaitescu & Meszaros 1999, who 
suggest that this phase should occur later). In addition, some 
steepening of the lightcurve decay may be due to a curved 
synchrotron spectrum: in fact at early times the spectral in- 
dex derived on the basis of BVRI photometric observations, 
de-reddened with Eb-v = 0.20, was flat {a = 0.61 ± 0.12, 
21.5 hours after the burst, Stanek et al. 1999b). Such a flat 
spectral index in the optical band must necessarily steepen 
at higher frequencies, to limit the emitted power. An esti- 
mate of such a steepening will come from the analysis of the 
X-ray afterglow flux, observed from 8 to 44.5 hours after 
the burst (Kuulkers et al. 1999). 

We conclude that the observed steepening of the 
lightcurve decay may then be the combined result of a 
curved synchrotron spectrum and of a coUimated fireball. 
In this scenario, the spectral index of the late (after 5 days) 
optical spectrum should be a ~ 1.1-1.2. 



4 DISCUSSION 

We have proposed a model to derive the amount of linear 
polarization observable from a coUimated fireball. We have 
shown that some degree of polarization can be observed even 
if the magnetic field is completely tangled in the plane of 
the fireball, as long as the fireball is observed slightly off- 
axis. One of the main virtues of the proposed model is in its 
easily testable predictions: i) the lightcurve of the degree of 
polarization has two maxima; ii) the observable polarization 
position angles are fixed between the first and the second 
maximum, being orthogonal between each other; iii) there 
is a strong link with the lightcurve of the total flux. 

The degree of polarization is predicted to be moderate, 
reaching 10% only if we are observing a coUimated fireball at 
its edge, and only for a short period of time. A larger degree 
of polarization would then suggest that the magnetic field is 
not completely tangled in the plane orthogonal to the line 
of sight, as suggested by Gruzinov & Waxman (1999). Up to 
now very few attempts have been done to measure linear po- 
larization in optical afterglows, and is therefore premature 
to draw any firm conclusion from the upper limit detected 
in GRB 990123 (Hjorth et al. 1999) and from the positive 
detection in the case of GRB 990510 (Covino et al. 1999, 
Wijers et al. 1999). Note however that the lightcurve be- 
havior of GRB 990510 matches our predictions, as well as 
the constant position angle of the observed polarization. At 
the time of the two first polarization measurements of GRB 
990510 (made two hours apart) the lightcurve of the total 
flux was decaying as t"^"' (Stanek et al. 1999a), i.e. it was 
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already steepening, in agreement with our model. Unfortu- 
nately, the polarization value measured one day later (Wijers 
et al. 1999) was not precise enough to further constrain the 
proposed scenario. 

It will be very interesting to explore in the future the as- 
sociation of a gradually steepening lightcurve and the pres- 
ence of polarization. We cannot exclude the possibility that 
some of the already observed afterglows can indeed be fltted 
by steepening power laws, but that the lack of data makes 
such an attempt meaningless. In addition, there are some op- 
tical afterglows (e.g. GRB 980326, Groot et al., 1998; GRB 
980519, Halpern et al., 1999) which showed a rapid decay. 
In these cases we may have observed only that part of the 
lightcurve corresponding to l/F > 9c + Oo- From these con- 
siderations we conclude that GRB 990510 may not be unique 
in its category and that a large fraction of gamma-ray burst 
afterglows can have some degree of optical linear polariza- 
tion. 
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